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Semiconductor  Measurement  Technology: 
Automatic  Determination  of  the  Interstitial  Oxygen  Content 
of  Silicon  Wafers  Polished  on  Both  Sides 

Warren  K.  Gladden,**  Stephen  R.  Slaughter, t 
Walter  M.  Duncan,^  and  Asian  Baghdadi* 

•National  Institute  of  Standards  and  Technology 
(formerly  National  Bureau  of  Standards) 
Gaithersburg,  MD  20899 

^Texas  Instrtunents,  Incorporated 
Dallas,  TX  75265 

ABSTRACT 

This  Special  Publication  contains  FORTRAN  and  PASCAL  computer  programs  which 
implement  sua.  ASTM  test  method  for  the  automatic  determination  of  the  interstitial  oxygen 
content  of  silicon.  The  programs  are  to  be  used  as  illustrative  examples  by  programmers 
wishing  to  implement  the  ASTM  algorithm  on  their  computers.  The  Publication  also 
includes  sample  data  that  can  be  used  to  test  the  computer  programs.  The  sample  data 
are  included  in  two  forms:   in  print,  and  on  an  MS-DOS  floppy  disk. 

Key  words:  algorithm;  analysis  of  infrared  spectra;  computer  programs;  interstitial  oxygen; 
silicon. 

I.  INTRODUCTION 

This  work  is  an  outgrowth  of  an  American  Society  for  Testing  and  Materials  (ASTM) 
project  to  automate  test  methods  for  the  determination  of  the  oxygen  and  carbon  con- 
tent in  silicon  using  IR  absorption.  The  ASTM  task  force  on  this  project  recommended 
a  test  procedure  for  obtaining  the  IR  spectra  and  an  algorithm  for  computing  the  oxygen 
content  from  the  spectra.  This  test  procedure  has  been  proposed  for  adoption  by  ASTM 
Committee  F-1  as  a  Standard  Test  Method.  This  Special  Publication  includes  as  examples 
a  FORTRAN  computer  program  and  a  PASCAL  program,  each  of  which  implements  the 
ASTM  aJgorithm.  The  FORTRAN  program  was  written  at  the  National  Institute  of  Stan- 
dards and  Technology,  and  the  PASCAL  progran  was  written  at  Texas  Instruments.  Also 
included  are  example  spectra  that  can  be  used  to  test  any  program  written  implementing 
the  ASTM  algorithm. 

*  Now  at  AT&T  Bell  Laboratories,  AUentown,  Pennsylvania. 
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The  proposed  ASTM  method  can  be  used  to  analyze  IR  spectra  obtained  on  either  dis- 
persive (DIR)  or  Fourier- transform  (FT-IR)  infrared  spectrophotometers.  It  can  be  used 
with  samples  ranging  in  thickness  from  0.3  mm  to  2.5  mm.  The  samples  must  be  polished 
on  both  sides.  The  lower  resistivity  limits  are  0.05  fl-cm  for  n-type  and  0.5  0-cm  for 
p-type  samples.  There  are  no  upper  resistivity  limits.  The  oxygen  content  of  float-zoned 
silicon  has  been  measured  with  resistivities  well  above  100  fi-cm.  The  range  of  oxygen 
concentrations  measured  by  this  method  is  1  X  10^^  to  2  x  10^®  atoms/cm~^. 

II.  OUTLINE  OF  THE  ALGORITHM 

The  interstitial  oxygen  content  of  silicon  can  be  determined  by  measuring  the  net  peak 
height  of  the  oxygen  vibrational  absorption  at  1107  cm~^.§  The  transmittance  through 
the  sample,  taking  into  account  the  additional  intensity  due  to  mvdtiple  reflections  between 
the  parallel  polished  surfaces  of  the  sample,  is  given 


(1  -  Rfe-'^ 


and  rearrangement  of  this  equation  yields  the  absorption  coefficient  as 


1,    r  -0.49  +  VO.2401  +  0.36T2  ) 
«  =   ^  /' 

where  the  silicon  reflectivity,  R,  has  been  set  equal  to  0.30,  its  value  over  the  carrier  con- 
centrations and  frequencies  considered  here.  Figure  1  shows  a  typical  silicon  transmittance 
spectrum  and  figure  2  shows  the  region  about  the  oxygen  peak  where  the  transmittance 
data  have  been  converted  to  absorption  coefficient  using  eq  (2).  To  determine  this  peak 
height,  a  baseline  is  drawn  between  the  two  minima,  one  on  either  side  of  the  peak.  The 
peak  height  is  determined  at  1107  cm'"-',  and  the  baseline  absorption  is  determined  by 
finding  the  intercept  of  the  baseHne  at  1107  cm"-'.  The  net  peak  height  is  then  the  differ- 
ence between  these  two  values  of  absorption  coefficient.  A  calibration  constant  has  been 
determined  relating  this  net  peak  height  to  the  quantity  of  interstitial  oxygen  in  parts  per 
million  (atomic)  using  absolute  techniques  such  as  charge  particle  and/or  photon  activation 


§  See  ASTM  Test  Method  F  121,  1987  Annual  Book  of  ASTM  Standards,  V.  10.05,  p. 
247. 

^  See  ASTM  Standard  Practices  F  120,  1987  Annual  Book  of  Standards,  V.  10.05,  p. 
245. 
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TRANSMITTANCE  SPECTRUM 
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Figure  1:  Transmittance  spectrum  of  a  p-type  Czochralski  silicon  wafer  with  a  resistivity 
of  9  n-cm  and  an  oxygen  content  of  20  parts  per  million  (atomic). 
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ABSORPTION  COEFFICIENT  SPECTRUM 
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Figure  2:  Absorption  coefficient  spectrum  of  a  p-type  Czochralski  silicon  wafer  with  a 
resistivity  of  9  17-cm  and  an  oxygen  content  of  20  parts  per  million  (atomic).  The  thick 
sections  of  the  line  indicate  the  wavenimiber  regions  used  to  calcvdate  the  oxygen  content 
(see  page  2). 
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analysis.  This  constant  is  6.28  cm-ppma."^ 

This  procedure  was  converted  to  a  computer  algorithm.  The  following  is  an  outline  of  the 
algorithm. 

•  The  measured  IR  transmittance  of  the  sample  is  read  from  a  data  file  and  only  those 
transmittances  between  1000  and  1300  wavenumbers  are  stored.  Three  regions  about 
the  interstitial  oxygen  minimum,  1010  to  1060  cm~^,  1200  to  1260  cm"-',  and  1090 
to  1123  cm~^,  are  considered  for  further  analysis.  The  first  two  regions  are  used  to 
determine  the  baseline  for  the  minimum.  By  the  method  of  least-squares,  a  third- 
order  polynomial  is  fitted  to  the  transmittance  data  in  each  of  these  two  regions.  For 
the  absorption  majcimum,  a  fourth-order  polynomial  is  used. 

•  The  two  transmittance  majcima,  the  minimum,  and  the  wavenumbers  at  which  they 
occur  are  determined  from  the  least-squares-fit  curves.  The  transmittance  values  are 
denoted  by  Tmaxi,  Tp^ak-,  and  Tmaxi-  The  respective  wavenumbers  are  denoted  by 
^  max  \i  CTpeak,  and 

^  max2  • 

•  The  baseline  transmittance  is  determined  from  a  linear  interpolation  between  the  two 
majcima  of  the  form 

T2  -  Ti 

Tba»e{o-)  =  Ti  -I  (o-  -  Cri). 

<T2  —  <Ti 

•  ocpeak  and  abase  are  determined  by  solving  eq  (2)  using  the  appropriate  calculated 
transmittance. 

•  The  net  absorption  coefficient  is  calculated  by 

+  A.  Baghdad!,  W.  M.  Bullis,  M.  C.  Croarkin,  Li  Yue-zhen,  R.  I.  Scace,  R.  W.  Series, 
P.  Stallhofer,  and  M.  Watanabe,  to  be  pubUshed  in  the  Journal  of  the  Electrochemical 
Society. 
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Lattice  absorption  at  1107  cm  ^  due  to  lattice  bands  contributes  an  additional  net  ab- 
sorption of  0.5  cm~^,^  so  the  oxygen  content  is  calculated  by 

ppm    atomic  =  6.28  x  {anet  —  0.5). 


This  algorithm  has  been  implemented  in  FORTRAN  and  PASCAL  computer  programs, 
which  are  listed  in  the  sections  III  A  and  IIIB,  respectively.  These  programs  should  only 
be  considered  as  examples  rather  than  as  the  approved  methods  for  implementing  the 
algorithm.  Section  IV  contains  four  sample  transmission  spectra.  These  spectra  (plus 
two  more)  have  also  been  written  on  5^-in.  MS-DOS*  floppy  disks  that  can  be  found 
in  a  pocket  on  the  end  of  this  Special  Publication.  Section  V  is  a  Hsting  of  the  output 
from  the  FORTRAN  and  PASCAL  programs  listed  in  section  III,  as  applied  to  all  six 
sample  spectra.  The  thicknesses  of  the  silicon  wafers  used  to  produce  the  IR  spectra 
are  included  in  the  output  listed  in  section  V.  The  oxygen  contents  determined  for  these 
spectra  should  match  the  oxygen  contents  listed  in  this  output.  Thus,  these  spectra  can 
be  used  to  test  a  computer  program  based  upon  the  algorithm  by  comparing  the  oxygen 
content  computed  using  the  program  to  the  listed  oxygen  content.  The  oxygen  contents 
in  section  V  are  listed  to  four  significant  figures,  rather  than  the  three  significant  figures 
which  is  appropriate  considering  the  precision  of  the  IR  measurement,  in  order  to  facilitate 
the  comparison  of  computer  programs. 


*  See  ASTM  Test  Method  F  121,  1988  Annual  Book  of  ASTM  Standards,  V.  10.05,  to 
be  published. 

*  Certain  commercial  equipment,  instruments  or  materials  are  identified  in  this  paper 
in  order  to  adequately  specify  the  experimental  procedure.  Such  identification  does  not 
imply  recommendation  or  endorsement  by  the  National  Institute  of  Standards  and  Tech- 
nology, nor  does  it  imply  that  the  materials  or  equipment  identified  are  necessarily  the 
best  available  for  the  purpose. 
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ANNOTATED  PROGRAMS 


Section  III  A.  FORTRAN  Program 


PROGRAM  ASTMDSP 

Q  Atik*******!******************  ASTMDSP 

C  *  * 

C  *  THIS  PROGRAM  COMPUTES  THE   INTERSTITIAL  OXYGEN  CONTENT  OF  * 

C  *  DOUBLE-SIDE  POLISHED  SILICON  SAMPLES   IN  PARTS  PER  MILLION  ATOMIC.  * 

C  *  THE  PROCEDURE   IS  A  MODIFICATION  OF  THAT  PRESENTED   IN  THE  GREY  * 

C  *  PAGES  OF  THE  ASTM  HANDBOOK.    INSTEAD  OF  TAKING  AVERAGE  VALUES  FOR  * 

C  *  THE  TWO  MAXIMA  AND  THE  MINIMUM,  A  LEAST-SQUARES  POLYNOMIAL   IS  FIT  * 

C  *  TO  SEVERAL  DATA  POINTS  AROUND  EACH  FEATURE.  ALSO,  THE  WAVENUMBERS  * 

C  *  WHERE  THE  MAXIMA  AND  MINIMUM  OCCUR  ARE  DETERMINED  FROM  THE  FIT.  A  * 

C  *  3RD-0RDER  POLYNOMIAL   IS  USED  '"OR  THE  MAXIMA  AND  A  4TH-0RDER  * 

C  *  POLYNOMIAL   IS  USED  FOR  THE  MINIMUM.  THE  4TH  ORDER  FIT  WAS  CHOSEN  * 

C  *   IN  ORDER  TO  BE  COMPATIBLE  WITH  THE  PROCEDURE  GIVEN   IN  THE  ANNEX  * 

C  *  FOR  THE  OXYGEN  DETERMINATION   IN  THIN  SLICES.  * 

C  *  * 

C  *  NOTE:  * 

C  *  IT   IS  COMMON  PRACTICE  TO  SPEAK  OF  THE   INTERSTITIAL  * 

C  *  OXYGEN  PEAK.  THIS   IS  ONLY  TRUE  WHEN  WORKING   IN  * 

C  *  ABSORBANCE.  ALL  DATA  REDUCTION   IN  THIS  PROGRAM   IS  * 

C  *  DONE   IN  TRANSMITTANCE,  THUS  WE  SPEAK  OF  THE  * 

C  *  INTERSTITIAL  OXYGEN  MINIMUM.  * 

C  *  * 

Q       ]|t4:]|C]fc>|e  lie  stole  sic  ******  4:  lie  sic  4c***  sec****    ASTMDSP  ********************************* 

DIMENSION  ABSRB(500) ,FUNC(200) , IFLAG(6) ,SIGMA(500) ,SPEC(600) ,W(6) , 
+    WAVE(600) ,WEIGHT(200) ,WORK(500) 
CHARACTER*40  F I LE I N , F I LEOUT 
DOUBLE  PRECISION  ANSWER(200) 
INTEGER  DAY, YEAR 

REAL  MAX1 ,MAX2,MINPK,MAXCAL,MINCAL,NET,NETABS,NODE(200) 
DATA  I  FLAG, W, WE  I GHT/6*0,1 010. ,1060. ,1090. ,1123. ,1200. ,1260. ,  . 
+  200*1./ 

1001  F0RMAT(3I2) 

1002  F0RMAT(A40) 

1003  F0RMAT(F10.4) 

1011  FORMATCIHI/'  ASTMDSP    ',12,'/', 12.'/'. 12, 

+  35X,'PAGE  ',11) 

1012  FORMATCIHI/'  ASTMDSP   ',12,'/', 12.'/', 12, 

+  34X,'PAGE  ',12) 

1013  FORMATC  *****spectrum*****  '  ,A40) 

1014  FORMATC  thickness  Is  '  ,f6.4, 'cm' ) 

1015  FORMATC/'  TC',F6.1,')  -  ',F6.4,'  residual  std.  dev.-',e10. 
+3/'  TC',F6.1,')  -  ',F6.4,'  residual  std.  dev.-' ,e10.3/'  T 
+C',F6.1,')  -  ',F6.4,'              residual  std.  dev.-' ,e10.3) 

1016  FORMATC/'     BASELINE  TRANSMITTANCE  AT   ',F6.1,'  WAVENUMBERS- ', F7 . 4) 

1017  FORMATC  — — >     ALPHACBASE)  -  ',F5.3,'         ALPHACPEAK)  -  ',F7.5/) 

1018  FORMATC     NET  ABSORBANCE  PEAK  HEIGHT  -  ',F8.6/) 

1019  FORMAT C 1 0X ,'*******************************'/, 1 0X ,' *  OXYGEN  CONTEN 
+T   IS   '   F5.2   '   PPMA  *'/  1 0X   '  ******sic************************ '  / ) 

NLINES-0 
IPAGE-1 

C 

Q  ********************************************************************* 

C    *     FILE   'MODEL. RES'    IS  USED  FOR  STORAGE  OF  THE  RESULTS.  THE  OUTPUT  * 

C    *     IS  PRINTED  TO  THIS  FILE   IN  SUCH  A  WAY  THAT  YOU  CAN  OBTAIN  A  * 

C    *     FORMATTED  HARD-COPY  OF  THE  FILE  AT  A  LATER  TIME.  * 
********************************************************************* 
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OPEN(UNIT-30,F ILE-'  MODEL. RES' , STATUS- ' NEW ) 

C 

Q  ****************************************************** 

C     *     INPUT  DATE,  NAME  OF  TRANSM I TTANCE  FILE,  AND  THICKNESS  OF  WAFER.  * 

c 

WRITE(6,*)    '     TODAY"S  DATE  ' 

READ (5,1 001 )MONTH, DAY, YEAR 
5  WRITE (6,*)    '     NAME  OF  DSP  TRANSM I TTANCE  DATA  FILE.' 
READ(5,1002)FILEIN 

WRITE(6,*)    '     THICKNESS  OF  WAFER.' 
READ(5,1003)TAU 

0PEN(UNIT-7,F ILE=F ILE IN,STATUS='OLD' ) 
IFCNLINES.EQ.O.AND. IPAGE.EQ.1)  THEN 
WR ITE(30,1011 )MONTH,DAY,YEAR, I  PAGE 
NL INES=NL INES+1 
END  IF 
WRITEOO,*)' 
WRITE (30,*)' 
WR ITE(30,1013)F ILE IN 
WR ITE(30,1014)TAU 

C 

C  *  READ  TRANSM I TTANCE  DATA  AND  SELECT  THOSE  DATA  POINTS  BETWEEN  1000  * 
C  *  AND  1300  WAVENUMBERS  (INCLUSIVE).  * 
C     *  FORMAT  OF  DATA  FILE:  * 


C  *  (1)  FIRST  LINE  CONTAINS  THE  NUMBER  OF  LINES  OF  DATA 

C  *  (2)  EACH  SUBSEQUENT  LINE  CONTAINS  THE  COUNTER  VALUE 

C  *  CORRESPONDING  TO  THE  LINE  OF  DATA,  WAVENUMBER ,  AND 

C  *  TRANSM I TTANCE. 


c 

K-1 

READ(7,*)NCNT 
DO  10   1-1 ,NCNT 

READ(7,*,ERR-115)LL,WAVENM, TRANSM 

I F ( ( WAVENM . GE . 1 000 . ) . AND . ( WAVENM . LE . 1 300 . ) )  THEN 
WAVE(K)-WAVENM 
SPEC(K)-TRANSM 
K-K+1 
GO  TO  10 
ELSE   IF  (WAVENM. GT.1 300. )  THEN 

GOTO  15 
END  IF 
10  CONTINUE 
15  NDP-K-1 

C 

C     *  DETERMINE  WAVENUMBER  REGIONS  FOR  LEAST-SQUARES  FIT  OF  * 

C     *         DATA  AT  MAXIMA  AND  MINIMUM  ABOUT  1107  WAVENUMBER  ABSORPTION  * 

C 
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MAXDEG-3 
KNT-0 

INDEX— 1 

DO  100  1-1 ,NDP 

IF(WAVE( I ) .GE.W(I))  INDEX-2 

IF(WAVE( I ) .GE.W(3))  INDEX-4 

IF(WAVE( I ) .GE.W(5))  INDEX-S 

IF( INDEX. LT.O)GOTO  100 

C 

C    *  CREATE  ARRAYS  OF  NODES  AND  FUNCTION  VALUES  FOR  LEAST  SQUARES  FIT.  * 

Q  3|C3|c)tcJ|(34c]K3lc3K4e3k3|e}|c3|e)k3|cilc3|c)Kik)lc)|e}|c>|c4c9|e4c>lc3k)lc4c)|c3k>^ 

C 

IF( IFLAGC INDEX) .EQ.1 )GOTO  100 
25   IF(KNT.EQ.O)  ITST-INDEX 
IF( ITST.NE. INDEX)  THEN 
KNT-0 
GO  TO  25 
END  IF 

35   IF(WAVE( I ) .LE.W( INDEX))  THEN 
KNT-KNT+1 
NODE(KNT)-WAVE( I ) 

IF( INDEX.EQ.4)N0DE(KNT)=N0DE(KNT)-1107. 
FUNC(KNT)-SPEC( I ) 
GO  TO  100 
END  IF 

C 

C    *  COMPUTE  LEAST  SQUARES  FIT.  * 

C 

45   IFLAGC I NDEX)-1 
EPS— 1  . 

IF( INDEX.EQ.4)  THEN 
MAXDEG-4 

CALL  LEAST(KNT , NODE , FUNC , WE  I GHT , EPS .MAXDEG , NDEG .WORK .ANSWER ) 
ll-INDEX-1 
GOTO  65 
ELSE 

MAXDEG=3 

CALL  LEAST ( KNT . NODE . FUNC .WE  I GHT , EPS .MAXDEG , NDEG , WORK . ANSWER ) 
ll-INDEX-l 

WAVCAL-QKPR IM(WORK.MAXDEG) 

IF(WAVCAL.GE.W( I  I ) . AND . WAVCAL . LE . W( INDEX))GOTO  65 
END  IF 

C 

C  *  ROUTINE  TO  USE  INSTEAD  OF  DERIVATIVE  PROCEDURE.  THIS  SECTION  OF  * 
C  *  CODE  IS  NECESSARY  SINCE  IT  IS  POSSIBLE  FOR  THE  FIT  NOT  TO  HAVE  A  * 
C    *     DISTINCT  Ml N  OR  MAX.  * 

c 
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STEP-(WAVE(2)-WAVE(1 ))/5. 

WAVNUM=W( I  I ) 

KK-0 

55   IF(WAVNUM.LE.W( INDEX))  THEN 
KK-KK+1 

ABSRB(KK)-EVAL(WAVNUM,MAXDEG,WORK,MAXDEG) 
SIGMA(KK)-WAVNUM 
WAVNUM-WAVNUM+ST  E  P 
GO  TO  55 
END  IF 

WAVCAL-MAXCAL ( S I GMA , ABSRB , KK ) 

C 

C    *  COMPUTE  RESIDUAL  STANDARD  DEVIATION  OF  THE  FIT  AND  THE  * 

C    *  WAVENUMBERS  AT  WHICH  MAXIMA  AND  MINIMUM  OCCUR.  * 

c 

65  SUM-0. 

DO  70  J-1 ,KNT 

SUM-SUM+ ( FUNC ( J )-SNGL ( ANSWER ( J ) ) ) **2 
70  CONTINUE 

RESID-SQRT(SUM/FL0AT(KNT-1 )) 

IF( INDEX. LT.1 .OR. INDEX. GT. 6)  THEN 

GOTO  135 
ELSE   IF( INDEX. EQ. 2)  THEN 
WAVE1-WAVCAL 

MAX 1 -EVAL ( WAVE 1 .MAXDEG .WORK .MAXDEG) 

R1-RESID 

GO  TO  100 
ELSE   IF( INDEX.EQ.4)  THEN 

CALL  PEAKCL (MAXDEG .WORK, WAVEP .M I NPK) 

R2-RESID 

GO  TO  100 
ELSE   IF( INDEX. EQ. 6)  THEN 

WAVE2-WAVCAL 

MAX2-EVAL(WAVE2 .MAXDEG .WORK .MAXDEG) 
R3-RESID 
ELSE 

GO  TO  135 
END  IF 
100  CONTINUE 

WRITE(30,1015)WAVE1 .MAX1 ,R1 .WAVEP .M I NPK . R2 , WAVE2 ,MAX2 , R3 

C 

C  *  CALCULATE  THE  BASE  TRANSM I TTANCE  AT  1107  WAVENUMBERS.  AND  THE  * 
C    *    ABSORPTION  COEFFICIENTS  AT  THE  PEAK  AND  BASE  OF  THE  ABSORPTION.  * 

c 

BASE- ( ( WAVE  2-WAVEP ) *MAX 1  +  ( WAVEP-WAVE 1 ) *MAX  2 ) / ( WAVE  2-WAVE 1 ) 
ALFAB-DSPCOR(BASE ,TAU) 
ALFAP-DSPCOR (M I NPK , TAU ) 
ALFNET-ALFAP-ALFAB 
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C    *  OUTPUT  REMAINING  RESULTS.  * 

c 

NETABS-ALFNET*TAU/2 . 303 
PPM-(ALFNET-0.5)*6.28 
WR I TE ( 30 , 1 01 6)WAVEP , BASE 
WR ITE(30,1017)ALFAB,ALFAP 
WRITE(30,1018)NETABS 
WRITE(30,1019)PPM 

C 

C    *  RESET   'I  FLAG'  ARRAY.  * 

c 

DO  110  MM-1 ,6 
110  IFLAG(MM)-0 

C 

Q  )|C}|C9k9|C3|C9fc9|e9|c4C3k3|C3k3|C3|C)|C4C3K9k3|e3|c:|(4C3lC3le3k4C3k3|C3fc)k3lC3k 

C    *  CHECK  FOR  FULL  PAGE.  * 

Q        3|e  4c    }k  )lc  )|C  9|C  3|C  3|c  lie  3|c  4c  ik  3k  i|c  A  )|c  3lc  3|C  3|C  3|e  i|c  lie  4c  lie  3lc  #  Ik  A  sfc  3k  4c  3|c  4c  ^ 

c 

NLINES-NLINES+16 

IF((NLINES+16) .GE.60)  THEN 
NLINES-1 
IPAGE-IPAGE+1 

IF( IPAGE.LT.10)  THEN 

WR I TE ( 30 , 1 01 1 )MONTH , DAY , YEAR , I  PAGE 
ELSE 

WR I TE ( 30 , 1 01 2 )MONTH , DAY , YEAR , I  PAGE 
END  IF 
END  IF 
GO  TO  5 

C 

Q  3k3k3k3k9k9k3k3k3k3k3k3k3k3k9k9k9k3k)k3k)k3k)kik3k3k3k3k3k3k3k3k3k3k3k3k3k3k3k3k3k3k3k3k3^ 

C    *  ERROR  MESSAGES  AND  TERMINATION.  * 

Q  )k)k9k)k)k3k)k3k)k3k3k3k3k)k3k>kik3k3k3k9k3k9k>k)k9k3k3kik3k3k3k3k)k)k3k3k3k}k)kik3k3k3k^ 

c 

115  WRITE(6,*)   '     ERROR  ENCOUNTERED   IN  READING  DATA.' 
GO  TO  200 

125  WRITE(6,*)   '     WAVENUMBER  CALCULATED   IS  OUTSIDE  OF  RANGE.' 

WR ITE(6,*) INDEX 

NUM-3*MAXDEG 

DO  130  JJ-1 ,NUM 
130  WRITE(6,*)JJ,W0RK(JJ) 

GO  TO  200 

135  WRITE(6,*)    '     HOW  DID  YOU  CHANGE  "INDEX"??????' 
200  STOP 
END 
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FUNCTION  DSPC0R(A1 ,A2) 

C 

C     *  THIS  FUNCTION  COMPUTES  ABSORPTION  COEFFICIENT  BASED  ON  * 

C    *  THE  DSP  TRANSMITTANCE  EQUATION.  * 

0 

TERM=(-0.49+SQRT(0.2401+0.36*A1**2))/(0.18*A1) 

DSPCOR— AL0G1  0(TERM)*2  .  303/A2 

RETURN 

END 

REAL  FUNCTION  MAXCAL ( ARRAY1 , ARRAY2 , NUM) 

C 

C     *  THIS  FUNCTION  FINDS  THE  MAXIMUM  TRANSMITTANCE  WITHIN  THE  * 

C     *  WAVENUMBER  RANGE  SPECIFIED  BY  THE  W( I  I )  AND  W( INDEX).  * 

Q       **********  *4:**«******«***«*************»»*««*****]|c*«ik*«**»«**«**]te**** 

0 

DIMENSION  ARRAY1 (*) ,ARRAY2(*) 
PEAK-0.25 
DO  10  KNT-1 ,NUM 
IF(ARRAY2(KNT) .LT.PEAK)GOTO  10 
MAXCAL-ARRAY1 (KNT) 
PEAK-ARRAY2(KNT) 
10  CONTINUE 
RETURN 
END 

REAL  FUNCTION  MINCAL(D) 

C 

C  *  WHEN  A  FOURTH-ORDER  POLYNOMIAL  IS  FIT  TO  THE  OXYGEN  TRANSMITTANCE  * 
C  *  MINIMUM,  THE  RESULTING  DERIVATIVE  IS  A  CUBIC  EQUATION.  THE  * 
C  *  THE  SOLUTION  OF  THIS  EQUATION,  IN  THE  SUBROUTINE  'PEAKCL',  IS  3  * 
C  *  ROOTS.  THIS  FUNCTION  FINDS  THE  ROOT  CLOSEST  IN  ABSOLUTE  VALUE  TO  * 
C  *  1107  WAVENUMBERS.  THIS  ROOT  SHOULD  BE  THE  WAVENUMBER  OF  THE  * 
C    *   INTERSTITIAL  OXYGEN  MINIMA.  * 

0 

DIMENSION  D(3) 

VALUE-ABS(D(1)) 

IK-1 

DO  20  KK-2,3 

I F(ABS(D(KK) ) .LT.VALUE)  IK-KK 
20  CONTINUE 

MINCAL-D( IK) 

RETURN 

END 
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SUBROUTINE  PEAKCL (NDEG .ARRAY ,WAVENM .VALUE) 

C 

C  *  CODE  TO  DETERMINE  COEFFICIENTS  OF  THE  FOURTH-ORDER  POLYNOMIAL  * 
C  *  FITTED  TO  THE  OXYGEN  MINIMUM  AND  THEN  SOLVE  THE  CUBIC  EQUATION  * 
C     *  FOR  THE  WAVENUMBER  AT  THE  MINIMUM.  * 

Q       :|:  4c  3lc  4c  3fc  3|c  3|c  3lc  3lc  3lotc  4!  >ic  sfc  4c  4: 3fc  4: 

C 

REAL  MINCAL 
DIMENSION  ARRAY(*) 
C0EFX3-4.*ARRAY(12) 

C 

SUM=0. 

DO  100  JJ=4,1 ,-1 
100      SUM=SUM+ARRAY( JJ) 

C0EFX2=-3 . *ARRAY( 1 2)*SUM+3 . *ARRAY( 1 1 ) 

C 

SUM1-0. 
SUM2-0. 
SUM3=0. 

DO  200  JJ=4,2,-1 

SUM2=SUM2+ARRAY(JJ+3) 

SUM3-SUM3+ARRAY( JJ-1 ) 

JJM1-JJ-1 

DO  200  LL=JJM1 ,1 ,-1 
200       SUM1-SUM1+ARRAY( JJ)*ARRAY(LL) 

COEFX 1  —2  .  *ARRAY (  1  2 ) * ( SUM1 -SUM2 )-2  .  *ARRAY (11)  *SUM3+2  .  *ARRAY (10) 

C 

SUM1-0. 
SUM2-0. 
SUM3-0. 
SUM4-0. 

DO  300  JJ-4,3,-1 

SUM2-SUM2+ARRAY(JJ-2) 

SUM3-SUM3+ARRAY( JJ+2) 

JJM1-JJ-1 

DO  300  LL»=JJM1  ,2,-1 
JJM2-JJ-2 

DO  300  MM=JJM2,1 ,-1 
300      SUM1=ARRAY( JJ)*ARRAY(LL)*ARRAY(MM) 
DO  400  LL-3,2,-1 
LLM1-LL-1 

DO  400  MM-LLM1 ,1 ,-1 
400      SUM4-SUM4+ARR AY ( LL ) *ARRAY (MM ) 

COEFXO— ARRAY(12)*(SUM1-(ARRAY(4)*(ARRAY(5)+ARRAY(6)) )+ 

*  ARRAY(3)*ARRAY(5)+ARRAY(2)*ARRAY(7)+ 

*  (ARRAYd )*(ARRAY(6)+ARRAY(7))))+ARRAY(11 )*(SUM4-SUM3)- 

*  ARRAY(10)*SUM2+ARRAY(9) 

C 

A2-C0EFX2/C0EFX3 
A1=C0EFX1/C0EFX3 
A0=COEFX0/COEFX3 


14 


CALL  CUBIC(A0,A1 ,A2,R1 ,R2,R3) 

R00T-MINCAL(R1 ,R2,R3) 

WAVENM-1107.+R0OT 

VALUE-EVAL ( ROOT , NDEG , ARRAY , NDEG ) 

RETURN 

END 


FUNCTION  QKPR IM(ARRAY,NUM) 

0 

C    *  THIS  FUNCTION  COMPUTES  THE  WAVENUMBER  AT  WHICH  THE  FIRST  * 

C    *  DERIVATIVE  OF  THE  3RD-0RDER  POLYNOMIAL  LEAST  SQUARES  FIT   IS  ZERO.  * 

Q       4c4c***************4c  ****************************  ***alc****:k********4c***» 

C 

DIMENSION  ARRAY(*) 

A-3.*ARRAY(9) 

SUM-O. 

DO  10  I l-NUM,1 ,-1 
10  SUM-SUM+ARRAY( I  I ) 

B— (2.*SUM*ARRAY(9)-2.*ARRAY(8)) 
SUM1-0. 

DO  30   I l-NUM,1 ,-1 
NUM1-I 1-1 

IF(NUM1 .NE.O)  THEN 
DO  20  KK-NUM1 ,1,-1 
20  SUM1-SUM1+ARRAY( I  I )*ARRAY(KK) 

END  IF 
30  CONTINUE 
SUM2-0. 
NUM2-2*NUM-2 
NUM1-NUM+1 
DO  40   I I-NUM1 .NUM2 
40  SUM2-SUM2+ARRAY( I  I ) 
SUM3-ARRAY(1 )+ARRAY(2) 

C-AR  R AY ( 9 ) * ( SUM1 -SUM2 ) -ARRAY ( 8 ) *SUM3+ARRAY ( 7 ) 
ARGMNT-(B**2)-4.*A*C 
IF(ARGMNT.GE.O)  THEN 
TERMI-SQRT(ARGMNT) 
QKPR IM-(-B+TERM1 )/(2.*A) 
IF(QKPR IM.LT.O. )QKPR IM-(-B-TERM1 )/(2.*A) 
RETURN 
ELSE 

QKPR  IM— 1  . 
RETURN 
END  IF 
END 
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SUBROUTINE  LEA.ST(M,X,F,W, EPS, MAXDEG,NDEG, ARRAY, R) 

C 

C    THE  SUBROUTINE  LEAST  AND  THE  FUNCTION  EVAL  CALCULATE  THE  LEAST 

C    SQUARES  POLYNOMIAL  APPROXIMATION  TO  A  SET  OF  DATA  SPECIFIED  BY  THE 

C    ARRAY  OF  M  NODES,  X,  WITH  CORRESPONDING  FUNCTION  VALUES  AND  WEIGHTS 

C    IN  THE  ARRAYS  F  AND  W,  RESPECTIVELY.  THE  WEIGHTS  MUST  ALL  BE 

C    POSITIVE.  THE  POLYNOMIAL  IS  DETERMINED  IN  LEAST  AND  EVALUATED  IN 

C    EVAL.  ON  INPUT  EPS  IS  THE  DESIRED  WEIGHTED  RMS  ERROR.  THE  CODE 

C    INCREASES  THE  DEGREE  OF  THE  FIT  IN  AN  ATTEMPT  TO  MEET  THIS  ERROR 

C    REQUEST.  ON  RETURN  EPS  IS  SF  TO  THE  WEIGHTED  RMS  ERROR  OF  THE  FIT. 

C    BECAUSE  EPS  IS  USED  FOR  BOTH  .NPUT  AND  OUTPUT,  IT  MUST  BE  A  VARIABLE 

C    IN  THE  CALLING  PROGRAM.  MAXDEG  IS  THE  HIGHEST  DEGREE  ALLOWED  AND  MUST 

C    BE  LESS  THAN  OR  EQUAL  TO  (M-1) .  THE  ACTUAL  DEGREE  OF  THE  FIT  IS 

C    RETURNED  IN  NDEG.  TO  FORCE  THE  CODE  TO  USE  THE  PARTICULAR  DEGREE 

C    MAXDEG,  SET  EPS  NEGATIVE  ON  INPUT.  THE  DOUBLE  PRECISION  VECTOR  R  GF  M 

C    WORDS  OUTPUTS  THE  DOUBLE  PRECISION  VALUES  OF  THE  POLYNOMIAL  FIT  AT 

C    EACH  OF  THE  DATA  POINTS  X(I).  THE  VECTOR  "ARRAY"  SPECIFIES  THE 

C    ORTHOGONAL  POLYNOMIAL  FIT  AND  PROVIDES  WORKING  STORAGE.  THE  DIMENSION 

C    OF  ARRAY  IN  THE  CALLING  PROGRAM  MUST  BE  AT  LEAST  2*M+3*MAXDEG .  THE 

C    ARRAYS  X,F,W,  ARRAY  AND  R  MUST  BE  DIMENSION  IN  THE  CALLING  PROGRAM. 

C 

DIMENSION  X(*)  ,F(*)  ,W(*)  ,ARRAY(0 
DOUBLE  PRECISION  R(*) , SUM, CK, TEMP 

C 

C    INITIALIZE  STORAGE  AND  CONSTANTS. 
C 

IB=MAXDEG+1 

IBL2=MAXDEG-1 

IC=IB+IBL2 

I0L1=IC+MAXDEG 

I1L1=I0L1+M 

RM=:M 

T0L=RM*EPS**2 

C 

C    CALCULATE  CONSTANT  FIT. 

C 

NDEG=0 
S=0. 

SUM,=0 .  ODO 
DO  1  1=1, M 
S=S+W(I) 

1  SUv.=SUW+DBLE(W(I))*D3LE(F(I}) 
RNO=S 

r 

C    CK  IS  THE  COEFFICIENT  C(0)  HERE. 
C 

CK=SUM/RNO 
ARRAY  (IC)=CK 
ERROR=0.0 
DO  2  1=1, M 
R(I)=CK 

2  ERR0R=ERR0R+W(I)*SNGL(CK-DB..E(F(I)))^^2 
IF  (NDEG. EQ. MAXDEG) GO  TO  14 
IF(EPS.LT.O.O)G0  TO  3 
IF(ERROR.LE.TOL)GO  TO  14 

C 

C    CALCULATE  LINEAR  FIT. 
C 

3  NDEG=1 
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ES=ERROR 
SUM=0 . ODO 
DO  4  1=1, M 

4  SUM=SUM+DBLE(W(I))*DBLE(X(I)) 

C 

C    CALCULATE  A(l) . 
C 

ARRAY (1)=SUM/RN0  ' 

C 

C    CALCULATE  Ql(.) . 
C 

S=0.0 
SUM=O.ODO 
DO  5  1=1, M 

ARRAY (IlLl+I) =X (I) -ARRAY  (1) 

S=S+W(I)*ARRAY(I1L1+I)*»2 

TEMP=DBLE(F(I))-R(I) 

5  SUM=SUM+DBLE(W(I)) *DBLE(ARRAY(I1L1+I)) *TEMP 
RN1=S 

C 

C    CK  IS  THE  COEFFICIENT  C(l)  HERE. 
C 

CK=SUM/RN1 
ARRAY (IC+1)=CK 

C 

C    CALCULATE  THE  VALUE  OF  THE  FIT  AT  THE  DATA  POINTS  AND 

C    ALSO  THE  RMS  ERROR. 

C 

ERROR=0.0 
DO  6  1=1, M 

R (I) =R (I) +CK*DBLE (ARRAY (IlLl+I) ) 

6  ERROR=ERROR+W (I) ♦SNGL (R (I) -DBLE (F  (I) ) ) *  *2 
IF(ERR0R.GT.ES.AND.EPS.GE.O.O)G0  TO  12 
IF(NDEG.EQ.MAXDEG)GO  TO  14 
IF(ERR0R.LE.T0L.AND.EPS.GE.O.O)G0  TO  14 
DO  7  1=1, M 

7  ARRAY (lOLl+I) =1.0 
NDEG=2 

K=2 

C 

C    GENERAL  FIT. 
C 

8  ES=ERROR 

C 

C    CALCUUTE  B(K) . 
C 

ARRAY (IBL2+K) =RN1/RN0 

C 

C    CALCULATE  A (K) . 
C 

SUM=O.ODO 
DO  9  1=1, M 

9  SUM=SUM+DBLE(W(I) ) ♦DBLE (X (I) ) *DBLE (ARRAY (IlLl+I)) **2 
ARRAY (K)=SUM/RN1 

C 

C    CALCULATE  QK(.)  OVERWRITING  ON  qK-2(.). 
C 

S=0.0 
SUM=O.ODO 
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DO  10  1=1, M 

ARRAY (lOLl+I) = (X (I) -ARRAY (K) ) *ARRAY (IlLl+I) 
1    -ARRAY (IBL2+K)*ARRAY(I0L1+I) 
S=S+W(I) *ARRAY(I0L1+I) **2 
TEMP=DBLE(F(I))-R(I) 

10  SUM=SLIM+DBLE (W (I)  )  ♦DBLE (ARRAY  (lOLl-^I) )  *TEMP 
RN0=RN1 

RN1=S 

C 

C    SWAP  INDICES  SO  II  REFERS  TO  STORAGE  OF  CIK(.) 

C    AND  10  TO  RK-1(.) . 

C 

IT=I0L1 

I0L1=I1L1 

I1L1=IT 

C 

C    CK  IS  THE  COEFFICIENT  C(K)  HERE. 
C 

CK=SUM/RN1 
ARRAY (IC+K)=CK 

C 

C    CALCULATE  THE  VALUE  OF  THE  FIT  AT  THE  DATA  POINTS  AND 

C    ALSO  THE  RMS  ERROR. 

C 

ERROR=0.0 
DO  11  1=1, M 

R  (I) =R (I) +CK*DBLE (ARRAY (IlLl+I) ) 

11  ERROR=ERROR+W (I) ♦SNGL (R (I) -DBLE (F (I) ) ) *  *2 
IF(ERR0R.GT.ES.AND.EPS.GE.O.O)G0  TO  12 
IF(NDEG.EQ.MAXDEG)GO  TO  14 
IF(ERR0R.LE.T0L.AND.EPS.GE.O.O)G0  TO  14 
NDEG=NDEG+1 

K=K+1 
GO  TO  8 

C 

C    HERE  IF  ERROR  INCREASED  ON  RAISING  DEGREE. 

C 

12  NDEG=NDEG-1 
ERROR=ES 

DO  13  1=1, M 

13  R  (I) =R (I) -CK*DBLE (ARRAY (IlLl+I) ) 

C 

C  EXIT. 
C 

14  EPS=SQRT(ERROR/RM) 
RETURN 

END 

FUNCTION  EVAL(Y,N, ARRAY, MAXDEG) 

C 

C    THE  FUNCTION  EVAL  EVALUATES  THE  ORTHOGONAL  POLYNOMIAL 
C    FIT  COMPUTED  BY  LEAST  AND  SPECIFIED  BY  THE  VECTOR 
C    ARRAY.  THE  FIT  OF  DEGREE  N  IS  EVALUATED  AT  THE  ARGU- 
C    MENT  Y.  N  MUST  BE  LESS  THAN  OR  EQUAL  TO  NDEG  AS 
C    RETURNED  FROM  LEAST.  LEAST  IS  CALLED  ONLY  ONCE  FOR 
C    EACH  FIT,  BUT  EVAL  IS  CALLED  ONCE  FOR  EACH  ARGUMENT 
C    AT  WHICH  WE  RQUIRE  THE  VALUE  OF  THE  FIT.  MAXDEG  MUST 
C    HAVE  THE  SAME  VALUE  AS  IN  THE  CALL  TO  LEAST. 
C 

DIMENSION  ARRAY (1) 
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IB=MAXDEG+1 
IC=MAXDEG+IB-1 

C 

C    EVALUATE  N=0,1  AS  SPECIAL  CASES. 
C 

IF(N.GT.O)G0  TO  1 

EVAL=ARRAY(IC) 

RETURN 

1  IF(N.GT.1)G0  TO  2 

EVAL=ARRAY  (IC)  +ARRAY  (IC+1)  *  (Y-ARRAY  (1) ) 
RETURN 

C 

C    GENERAL  RECURRENCE  RELATION. 
C 

2  DKP2=ARRAY(IC+N) 

DKP1=ARRAY (IC+N-1) ♦ (Y-ARRAY (N) ) *DKP2 
NL2=N-2 

IF(NL2.LT.1)G0  TO  4 
DO  3  L=1,NL2 
K=1+NL2-L 

DK=ARRAY (IC+K) + (Y-ARRAY (K+1) ) ♦DKPl 
1      -ARRAY  (IB-»-K)*DKP2 
DKP2=DKP1 

3  DKP1=DK 

4  EVAL=ARRAY (IC) + (Y-ARRAY (1) ) *DKP1 
1      -ARRAY (IB) *DKP2 

RETURN 
END 
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SUBROUTINE  CUBIC (AO, A1,A2,R1,R2,R3) 

C 

C    *  THIS  SUBROUTINE  WAS  WRITTEN  BY  JERRY  LOWNEY  * 

C    ♦  TO  SOLVE  THE  CUBIC  EQUATION:  * 

C    *  X*3  +  A2*X^2  *  A1*X  +  AO     =  0  * 

C     ^tf^f^itifi^t^ift^firm**********  ************  ************** 

c 

COMPLEX  X1,X2,X3,T,U,D,E,F 

P=A2 

Q=A1 

R=AO 

A=1./3*(3*Q-P**2) 

B=l . /27* (2*P**3-9*P*Q+27*R) 

D=B**2/4+A**3/27 

E=-B/2+CSQRT(D) 

IF(E.Eq.(0,0))  THEN 
T=0 

GO  TO  15 
END  IF 

IF(AIMAG(E) .EQ.O. .AND.REAL(E) .LT.O.)THEN 

T=-EXP (ALOG (CABS (E)  )  /3) 

GO  TO  15 
END  IF 
T=CEXP(CL0G(E)/3) 

15  F=-B/2-CSQRT(D) 

IF(F.EQ.(0,0))  THEN 
U=0 

GO  TO  16 
END  IF 

IF(AIMAG(F)  .EQ.O.  .AND.REAL(F)  .LT.O.)THEN 

U=-EXP (ALOG (CABS (F) ) /3) 

GO  TO  16 
END  IF 
U=CEXP(CL0G(F)/3) 

16  X1=:T+U 

X2=- . 5* (T+U) ♦ . 5* (T-U) ♦CSQRT ( (-3 , 0) ) 
X3=- . 5* (T+U) - . 5* (T-U) *CSQRT ( (-3 , 0) ) 

Rl=:Xl-P/3 

R2=X2-P/3 
R3=X3-P/3 
RETURN 
END 
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Section  IIIB.  PASCAL  Program 
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ASTMDSP     [Pascal  Version] 
program  as tmdsp ( Input , output , f n) ; 
(*F-1  Proposal  p  112*) 

(*Interstitial  Oxygen  Content  Of  Silicon  Slices  By*) 
(*Computer-Assisted  Infrated  Spectrophotometry*) 
(*10/16/85*) 

(*Double  Side  Polished*) 

(*  

*)   (*   Modified  by  Warren  Gladden 

(*   date:  8/21/87 

(*  

*) 


const 

maxr=32 ; 
maxc=9 ; 


(*max  data  points  in  fit*) 
(*max  number  of  variables  in  fit*) 


type 

sstring 

drayr 

drayc 

ddrayl 

ddray2 


packed  array [1.. 40]  of  char; 
array [ 1 . .maxr]  of  real; 
array [ 1 .. maxc ]  of  real; 
array [ 1 . .maxr , 1 . .maxc ]  of  real; 
array [ 1 . .maxc , 1 . .maxc]  of  real; 


(*filename  type*) 
(*array  types  for  fits*) 
(*  *) 
(*  *) 
(*  *) 


yar 


coef ,nal 

: drayc ; 

(*ls  fit  vars*) 

x.y.yemp, 

dual 

: drayr ; 

(*ls  vars*) 

xl : array [ 1 .. 1001 ]  of  real; 

(*wavenumber  array*) 

yl : array [1. .1001]  of  real; 

(*transmittance  array*) 

corco ,yco ,xco 

:  real 

(*ls  vars*) 

xint ,xwave 

:  real 

(*interpolate  vars*) 

smax2 , speak, tmaxl , smaxl , tmax2 

:  real 

(*astra  vars*) 

tr , thk 

:  real 

(*temp  trans ,  thick*) 

rsdmaxl , rsdmax2 , rsdpeak 

:  real 

(*residual  vars*) 

nmaxl , nmax2 , npeak 

:  real 

(*#  pts  fit*) 

apeak, abase , tpeak, temp , tempi 

:  real 

(*abs  pk.base*) 

ppm, anet , latnet 

:  real 

(*conc  var*) 

tbase, tbl,tb2 

:  real 

(*baseline  var*) 

al , a2 , a3 , alphas 

:  real 

(*abs  coef  var*) 

e , del , f irstx , rof f set 

:  real 

(*uncertainty  var*) 

sumy2 , ysum2 

:  real 

(*stan.  dev.  var*) 

nrow,ncol ,num,noffset 

: integer ; 

(*array  sizes  and  degree*) 

startwave , endwave , wavel , wave 2 

: integer ; 

(*wavenumber  limit  for  fit*) 

dummy , i , j , xtemp 

: integer ; 

(*loop  and  temp,  var*) 

f name , f names , j  unk 

: string [ 14 

] ; (*f ilename ,  ref,  samp*) 

fn 

: text; 

(*file  var*) 

(*  

*)  procedure  base; 
begin 

tbl :=smax2- speak; 


(*baseline  @  oxygen*) 
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tb2 : =speak- smaxl ; 
tbase :=tbl*tmaxl+tb2*tmax2 ; 
tbase :=tbase/(smax2-smaxl) ; 
end; 

(*  

*)  procedure  alpha;  (*absorption  coefficient*) 

begin 

al :=-0.49+sqrt(0. 2401+0. 36*sqr ( tr) ) ; 
a2:=al/(0.18*tr) ; 
a3:=ln(a2); 
alphas :=-l/thk*a3; 
end; 

(*  

*)  procedure  uncertainty;  (*non-s3mimetrical  uncertainty*) 

begin 

e:  =  (rsdmaxl/(ninaxl-4)  )  +  (rsdpeak/(npeak-5)  )  +  (rsdniax2/(ninax2-4) ) ; 
e:=2*sqrt(e) ; 
end; 

(*  

*)  procedure  standev;  (*standard  deviation*) 

begin 

svuny2 :  =0 ; 

ysum2 :=0 ; 

j:=0; 

for  i:=l  to  num  do 
begin 

if  (xl[i]>=2000)  and  (xl [ i]<=2060)  then 
begin 

sumy2 : =yl [ i ] *y 1 [ i ] +sumy2 ; 
ysum2 : =yl [ i ] +ysum2 ; 

j:=j+i: 

end ; 
end; 

ysum2  :=sqr (ysuin2) ; 
ysuin2  :=ysuni2/j  ; 
suiny2  :=s\iiny2-ysuin2 ; 
sumy2  :=suiny2/(  j  -2) ; 
sumy2 :=sqrt(sumy2) ; 
end; 

(*  

*)  procedure  OpenlnFile; 
var 

OK  : boolean; 

begin 
repeat 

write ('Enter  name  of  transmittance  file:  '); 

readln(fnaine)  ; 

assign(fn, fname) ; 

{$!-}  reset(fn)  {$1-); 

OK:=(IOresult  =  0); 

if  not  OK  then 

wr ite In (' Cannot  find  file:   ', fname); 
until  OK; 
end; 
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(*-   

*)  procedure  OpenFile; 
begin 

assign(fn,fnaine)  ; 

rewrite(fn) ; 
end; 

(*  

*)  procedure  sorttrans;  (*move  target  trans  values*) 

begin 

for  i:=l  to  num  do 

begin 

if  (xl [i]>=wavel)  and  (xl [ i]<=wave2)  then 
begin 

if  j=l  then  startwave :=i ; 
x[j]:=xl[i]; 

y[j]:=yi[i]; 
j:=j+i; 

end; 
end; 

nrow : =j - 1 ; 

endwave : =s tartwave+j - 2 ; 
end; 

(*  

*)  procedure  minmax;  (*sort  trans  values*) 

begin 

for  i:=l  to  nrow-1  do 
for  j :=i+l  to  nrow  do 
begin 

if  y[i]<y[j]  then 
begin 

temp:=y[i];  tempi :=x[i]; 
y[i]:=y[j];  x[i]:=x[j]; 
y[j]:=temp;  x[ j ] :=templ ; 
end; 
end; 

end; 

(*  

*)  procedure  dataout;  (*output  fit  results  to  file*) 


var 
i 


: integer ; 
to  ' ,wave2) ; 


begin 

writeln(fn,wavel , ' 
writeln(fn) ; 

writeln(fn, '  y  ycalc  residual'); 

writeln(fn) ; 

for  i:=l  to  nrow  do 

writeln(fn,x[i] :7:4, '  ' ,y [ i] : 10 : 6 , '  ' ,yemp [ i] : 10 : 6 , dual [ i] :9:3) ; 
writeln(fn) ; 

writeln(fn, ' correlation  coef  is  ' , corco : 8 : 5) ; 

writeln(fn) ; 

for  i:=l  to  ncol  do 

writeln(fn, ' coefficient     ' , coef [ i ] ) ; 
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writeln(fn) ; 
end; 

(*  

*)  procedure  s2 

(var  x:ddrayl; 
var  y:drayr; 
var  a:ddray2; 
var  g:drayc; 
var  nrow.ncol : integer) ; 


(*least  squares  fit  to  #  of*) 
(*x,y  points  in  nrow*) 


var 
i.k,l 


: integer; 


begin 

for  k:=l  to  ncol  do 
begin 

for  1:=1  to  k  do 
begin 

a[k.l]:=0; 
for  i:=l  to  nrow  do 
begin 

a[k,l]:=a[k,l]+x[i.l]*x[i.k]; 
if  kO>l  then  a[  1  ,k]  :=a[k,  1] 
end 

end; 
g[k]:=0; 

for  i:=l  to  nrow  do 
g[k]:=g[k]+y[i]*x[i,k] 

end 

end; 

(*  


*)  procedure  curvel  (var  b  :ddray2; 

y  : drayc ; 

var  coef  : drayc; 

ncol  : integer) ; 


var 
w 

cf  ile 

i,j ,k,l,nv,rw,cl,n,ll 

sure , triad , keep , sum , t , ab , great 
*  


: array [ 1 . .maxc , 1 . .maxc ]  of  real; 
; array [ 1 . .maxc , 1 . .maxc]  of  integer; 
: integer; 
; real ; 


( 

*)  procedure  xchng(var  a,b:real); 


var 
keep 

begin 

keep :=a; 

a:=b; 

b : =keep ; 
end; 

(*  

*)  procedure  curve2; 


: real ; 


(*swap*) 
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var 

i,j,k,l,ll  :integer; 
(*  

*)  procedure  curve3; 


var 

1  : integer ; 

begin 

if  rwOcl  then 
begin 

sure :=- sure ; 
for  1:=1  to  n  do 

xchng(b[rw,l] ,b[cl,l]) ; 
if  nv>0  then 

for  1:=1  to  nv  do 

xchng(w[rw, 1] ,w[cl , 1] ) 

end 

end; 

begin 
nv:=l ; 
n:=ncol ; 

for  i:=l  to  n  do 
begin 

w[i,l]:=y[i]; 
cfile[i,3] :=0 
end; 
sure:=l .0; 
for  i:=l  to  n  do 
begin 

great: =0,0; 
for  j :=1  to  n  do 
begin 

if  cfile[j  ,3]Ol  then 
begin 

for  k:=l  to  n  do 
begin 

if  cfile[k,3]<l  then 

if  abs(b[ j ,k] )>great  then 
begin 
rw:=j  ; 
cl:=k; 

great :=abs(b[j ,k] ) 
end 

end 

end 

end; 

cfile[cl,3] :=cfile[cl,3]+l; 
cf ile [ i , 1 ] :=rw; 
cfile[i,2] :=cl; 
curve 3 ; 

triad:=b[cl,cl] ; 
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sure : =sure*tr iad ; 
b[cl,cl] :=1.0; 


for  1:=1  to  n  do 

b[cl,l] :=b[cl,l]/triad; 
if  nv>0  then 

for  1:=1  to  nv  do 

w[cl,l] :=w[cl,l]/triad; 

for  11 :=1  to  n  do 
begin 

if  llOcl  then 
begin 

t:=b[ll,cl]; 
b[ll,cl] :=0.0; 
for  1:=1  to  n  do 

b[ll,l]:=b[ll,l]-b[cl,l]*t; 
if  nv>0  then 

for  1:=1  to  nv  do 

w[ll,l]:=w[ll,l]-w[cl.l]*t; 

end 

end 
end; 
end; 

begin 
curve 2 ; 

for  i:=l  to  n  do 
begin 

l:=n-i+l; 

if  cfile[l,l]Ocfile[l,2]  then 
begin 

rw:=cf ile [1,1] ; 
cl:=cfile[l,2] ; 
for  k:=l  to  n  do 

xchng (b[k,rw] ,b[k,cl]) 

end 

end; 

for  i:=l  to  n  do 
coef [i] :=w[i,l] ; 

end; 

(*  

*)  procedure  curve (x,y  :drayr; 

var  yemp  : drayr ; 

var  dual  : drayr ; 

var  coef  :drayc; 

var  nal  rdrayc; 

nrow  : integer; 

var  ncol  : integer); 

var 

xmatr  : ddrayl ; 

a  : ddray2 ; 

g  : drayc ; 
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i.j.nm  : integer; 

xi  ,yi  ,yc ,  srs ,  see ,  suni__y ,  sum_y2  :real; 

begin 

for  i:=l  to  nrow  do 
begin 

xi:=x[i]; 
xmatr[i, 1] :=1 .0; 
for  j :=2  to  ncol  do 

xmatr[i,  j  ]  :=xinatr[i,  j  -  l]*xi 
end; 

s2(xmatr ,y ,a,g,nrow,ncol) ; 
curvel(a,g,coef ,ncol) ; 
stim_y :  =0 . 0 ; 
sum_y2 : =0 . 0 ; 
srs:=0.0; 

for  i:=l  to  nrow  do 
begin 

yi:=y[i]; 

yc:=0.0; 

for  j :=1  to  ncol  do 

yc :=yc+coef [ j ]*xmatr [ i , j ] ; 
yemp [ i ] : =yc ; 
dual[i] :=yc-yi; 
srs :=srs+sqr (dual [ i] ) ; 
sum_y :  =suin_y+y  i ; 
siam__y2  :  =sum_y2+y i*y i ; 
end; 

corco:=sqrt(l  .0- srs/ (sum_y2-sqr(suin_y) /nrow) ) ; 
if  nrow=ncol  then  nm:=l 
else  nm:=nrow-ncol ; 
see :=sqrt(srs/nm) ; 
for  i:=l  to  ncol  do 

nal[i] :=see*sqrt(a[i, i] ) 
end; 

(*  

*)  procedure  polythree; 
begin 

xco :=x[ 1 ] -2 ; 

temp : =0 . 0 ; 

tempi :=0 . 0 ; 

for  i:=l  to  400  do 

begin 

xco :=xco+0 , 01 ; 

yco:=coef [1 ]+coef [2]*xco+coef [ 3]*sqr(xco)+coef [4]*xco*sqr(xco) ; 

if  yco>temp  then 

begin 

temp : =yco ; 
tempi :=xco; 

end;  (*if*) 
end;  (*loop*) 
xtemp : =trunc ( temp 1 ) ; 
xint : = temp 1 - trunc ( temp 1 ) ; 
xwave:=xl [xtemp+startwave]+xint*l .925; 
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end;   (*poly three*) 

(*  

*)  procedure  polyfour; 
begin 

xco:=x[nrow] -2 ; 

temp: =10.0; 

tempi : =0.0; 

for  i:=l  to  400  do 

begin 

xco:=xco+0.01 ; 

yco :=coef [1 ]+coef [2]*xco+coef [3]*sqr(xco)+coef [4]*xco*sqr(xco)+coef [5]*sqr(xco)* 

if  yco<temp  then 

begin 

temp : =yco ; 
tempi :=xco ; 
end;  (*if*) 
end;  (*loop*) 
xtemp : =trunc ( tempi ) ; 
xint : = temp 1 - t rune ( temp 1 ) ; 

xwave :=xl [xtemp+startwave]+xint*l . 925 ;  ' 
end;  (*polyfour*) 

(*  

*)  procedure  Isl; 
begin 

j:=i: 

wavel:=1010; 
wave2:=1060; 
sorttrans ; 

for  i:=l  to  nrow  do 
x[i]:=i; 
ncol:=3; 
begin 

ncol :=ncol+l ; 

curve (x , y , yemp , dual , coef , nal , nrow , ncol) ; 

dataout ; 
end; 
minmax ; 
polythree; 
smaxl :=xwave ; 
tmaxl :=temp ; 
rsdmaxl : =nal [ 1 ] ; 
nmaxl :=nrow; 
writeln(fn) ; 

writeln(fn, ' tmaxl  ', tmaxl); 

writeln(fn, ' smaxl  ', smaxl); 

writeln(fn, ' rsdmaxl  ', rsdmaxl); 

writeln(fn, '#  pts,  fit  ',nrow); 
writeln(fn) ; 
end;  (*lsl*) 

(*  

*)  procedure  ls2; 
begin 

j:=i: 

wavel  -.=1090; 
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wave2:=1123; 
sorttrans ; 

for  i:=l  to  17  do 
x[i]:=i; 
ncol :=4; 
begin 

ncol :=ncol+l ; 

curve (x,y,yemp, dual ,coef ,nal,nrow,ncol) ; 

dataout ; 
end; 
minmax ; 
poly four; 
speak :=xwave ; 
tpeak:=temp ; 
rsdpeak : =nal [ 1 ] ; 
npeak:=nrow; 
writeln(fn) ; 

writeln(fn, ' tpeak  '.tpeak); 
writeln(fn, ' speak  '.speak); 
writeln(fn, 'rsdpeak  ' , rsdpeak) ; 

writeln(fn, '#  pts.  fit  ',nrow); 
writeln(fn) ; 
end;  (*ls2*) 



*)  procedure  ls3; 
begin 

j:=i; 

wavel :=1200; 
wave2:=1260; 
sorttrans ; 

for  i:=l  to  32  do 
x[i]:=i; 
ncol : =3 ; 
begin 

ncol :=ncol+l ; 

curve (x.y.yemp, dual, coef,nal,nrow, ncol) ; 

dataout ; 
end; 
minmax ; 
polythree ; 
smax2 : =xwave ; 
tmax2 : =temp ; 
rsdmax2 : =nal [ 1 ] ; 
nmax2 :=nrow; 
writeln(fn) ; 

writeln(fn, ' tmax2  ',tmax2); 
writeln(fn, ' smax2  ',smax2); 
writeln(fn, ' rsdmax2  ',rsdmax2); 
writeln(fn, '#  pts.  fit  ',nrow); 
writeln(fn) ; 
end;  (*ls3*) 

(*  

*)  procedure  fileout;  (*output  results*) 

begin 
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fname:='resltd51 .dat' ; 
OpenFile ; 

writeln(fn, ' sample  file  name 
vn:iteln(fn,  '  thickness  (mils) 
writeln(fn) ; 
writeln(fn, 'tbase 
writeln(fn, 'tpeak 
writeln(fn, 'abase 
writeln(fn, 'apeak 
writeln(fn) ; 
writeln(fn, 'anet 
writeln(fn, 'uncertainty 
wr ite In ( f n, ' standard  deviation 
writeln(fn, 'ppm  (6.28) 
close(fn) ; 


, f names) ; 
,thk/0. 00254); 


, tbase) ; 
, tpeak) ; 
, abase) ; 
, apeak) ; 


, anet) ; 
,e) ; 

,sumy2) ; 
,ppm) ; 


end; 

(*-- 


*)  begin  (*MAIN  PROGRAM*) 
OpenlnFile ; 
f names : =fname ; 

readln(fn,num) ; 

writeln(num) ; 

for  i:=l  to  num  do 

begin 

readln(fn, dummy ,xl [ i] ,yl [ i] ) ; 

if  yl[i]>1.0  then  yl [i] :=yl [i]/100; 
end; 
close(fn) ; 
del:=xl[2]-xl[l] ; 
f irstx:=xl [ 1 ] ; 

roffset :=(f irstx/del) ;  (*  nof fset :=int (roff set)  ; *) 

writeln(rof f set , del , f irstx) ; 
writeln( ' thickness  (cm)  '); 
readln(thk) ; 
thk:=thk*1.0; 

f name :  = ' leas  t . da t ' ; 
OpenFile ; 

writeln(fn, ' sample  file  name  '.fnames); 
writeln(fn, ' thickness  (cm)  ',thk/1.0); 
writeln(fn) ; 


close(fn) ; 
base ; 

tr :=tbase; 
alpha; 

abase : =alphas ; 
tr :=tpeak; 
alpha ; 

apeak : =alphas ; 
anet : =apeak-abase ; 
latnet:=anet-0. 5 ; 


Isl; 
ls2; 
ls3; 
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ppm : =latnet*6 . 28 ; 

uncertainty; 

standev; 

writeln( 'OXYGEN  PPM    (6.28)  '.ppm); 
writeln; 

write In ('UNCERTAINTY  ',e); 
writeln; 

writeln(  '  STANDARD  DEVIATION     '  ,  sximy2)  ; 
f ileout ; 
end. 
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Section  IV.  SAMPLE  DATA 


The  data  are  listed  in  three  columns,  with  two  sets  of  three  columns  per  page.  The  first 
column  is  just  an  index  number.  The  second  column  lists  the  wavenumber,  and  the  third 
column  lists  the  transmittance  at  that  wavenumber.  The  single  integer  number  at  the  head 
of  each  spectrum  is  the  number  of  data  points  in  that  spectrum. 
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Samp  1 

le:  1 
840 

D07 

261 

501 

.401 

0 

.46960 

262 

503 

.330 

0 

.46700 

263 

505 

.258 

0 

.46442 

264 

507 

.187 

0 

. 46066 

265 

509 

.115 

0 

.45473 

266 

511 

.044 

0 

.44725 

267 

512 

.972 

0 

.44130 

268 

514 

.901 

0 

.44220 

269 

516 

.829 

0 

.44576 

270 

518 

.758 

0 

.44822 

271 

520 

.686 

0 

.45390 

272 

522 

.615 

0 

.45969 

273 

524 

.543 

0 

.46176 

274 

526 

.471 

0 

.46235 

275 

528 

.400 

0 

.46203 

276 

530 

.328 

0 

.46114 

277 

532 

.257 

0 

.46034 

278 

534 

.185 

0 

.45943 

279 

536 

.114 

0 

.45866 

280 

538 

.042 

0 

.45804 

281 

539 

.971 

0 

.45698 

282 

541 

.899 

0 

.45474 

283 

543 

.828 

0 

.45230 

284 

545 

.756 

0 

.45063 

285 

547 

.685 

0 

.44929 

286 

549 

.613 

0 

.44864 

287 

551 

.542 

0 

.44810 

288 

553 

.470 

0 

.44630 

289 

555 

.398 

0 

.44405 

290 

557 

.327 

0, 

. 44202 

291 

559 

.255 

0, 

.44026 

292 

561 

.184 

0, 

.43863 

293 

563, 

.112 

0, 

.43723 

294 

565, 

.041 

0, 

,43681 

295 

566, 

,969 

0, 

,43722 

296 

568, 

,898 

0. 

,43896 

297 

570, 

,826 

0. 

,44246 

298 

572. 

.755 

0. 

.44637 

299 

574. 

,683 

0. 

,44892 

300 

576. 

,612 

0. 

44954 

301 

578. 

,540 

0. 

44934 

302 

580. 

,469 

0. 

44890 

303 

582. 

,397 

0. 

44707 

304 

584, 

326 

0. 

44355 

305 

586. 

254 

0. 

43886 

306 

588. 

182 

0. 

43340 

307 

590. 

111 

0. 

42661 

308 

592. 

039 

0. 

41782 

309 

593. 

968 

0. 

40774 

310 

595. 

896 

0. 

39597 

311 

597. 

825 

0. 

38115 

312 

599. 

753 

0. 

36256 

313 

601. 

682 

0. 

34260 

314 

603. 

610 

0. 

32663 

315 

605. 

539 

0. 

31448 

316 

607. 

467 

0. 

30306 

317 

609. 

396 

0. 

29526 

318 

611 

.324 

0 

. 29407 

319 

613 

.253 

0 

.29794 

320 

615 

.181 

0 

.30223 

321 

617 

.109 

0 

.30424 

322 

619 

.038 

0 

.30586 

323 

620 

.966 

0 

.31074 

324 

622 

.895 

0 

.32582 

325 

624 

.823 

0 

.35559 

326 

626 

.752 

0 

.39533 

327 

628 

.680 

0 

.43551 

328 

630 

.609 

0 

.46502 

329 

632 

.537 

0 

.47996 

330 

634 

.466 

0 

.48596 

331 

636 

.394 

0 

.48837 

332 

638 

.323 

0 

.49016 

333 

640 

.251 

0 

.49234 

334 

642 

.180 

0 

.49368 

335 

644 

.108 

0 

.49471 

336 

646 

.036 

0 

.49555 

337 

647 

.965 

0 

.49591 

338 

649, 

.893 

0 

.49675 

339 

651, 

.822 

0, 

.49760 

340 

653, 

.750 

0, 

.49810 

341 

655, 

,679 

0, 

.49854 

342 

657, 

,607 

0. 

,49898 

343 

659, 

,536 

0, 

,49959 

344 

661. 

,464 

0. 

,49968 

345 

663. 

,393 

0, 

,49990 

346 

665. 

,321 

0. 

,49975 

347 

667. 

,250 

0. 

,49796 

348 

669. 

178 

0. 

49785 

349 

671. 

,107 

0. 

,50024 

350 

673. 

035 

0. 

.50163 

351 

674. 

963 

0. 

50216 

352 

676. 

892 

0. 

50297 

353 

678. 

820 

0. 

50352 

354 

680. 

749 

0. 

.50305 

355 

682. 

677 

0. 

50219 

356 

684. 

606 

0. 

50201 

357 

686. 

534 

0. 

50235 

358 

688. 

463 

0. 

50276 

359 

690. 

391 

0. 

50286 

360 

692. 

320 

0. 

50254 

361 

694. 

248 

0. 

50323 

362 

696. 

177 

0. 

50546 

363 

698. 

105 

0. 

50767 

364 

700. 

034 

0. 

50817 

365 

701. 

962 

0. 

50730 

366 

703. 

890 

0. 

50631 

367 

705. 

819 

0. 

50468 

368 

707. 

747 

0. 

50161 

369 

709. 

676 

0. 

49768 

370 

711. 

604 

0. 

49243 

371 

713. 

533 

0. 

48546 

372 

715. 

461 

0. 

47828 

373 

717. 

390 

0. 

47242 

374 

719.318 

0. 

46809 

375 

721. 

247 

0. 

46433 

376 

723. 

175 

0. 

46016 

34 


377 

725 

.104 

0 

.45587 

378 

727 

.032 

0 

.45224 

379 

728 

.961 

0 

.44886 

380 

730 

.889 

0 

.44554 

381 

732 

.817 

0 

.44276 

382 

734 

.746 

0 

.43980 

383 

736 

.674 

0 

.43686 

384 

738 

.603 

0 

.43605 

385 

740 

.531 

0 

.43816 

386 

742 

.460 

0 

.44186 

387 

744 

.388 

0 

.44549 

388 

746 

.317 

0 

.44848 

389 

748 

.245 

0 

.45117 

390 

750 

.174 

0 

.45307 

391 

752 

.102 

0 

.45431 

392 

754 

.031 

0 

.45566 

393 

755 

.959 

0 

.45651 

394 

757 

.888 

0 

.45704 

395 

759 

.816 

0 

.45794 

396 

761 

.745 

0 

.45926 

397 

763 

.673 

0 

.46062 

398 

765 

.601 

0 

.46116 

399 

767 

.530 

0 

.46109 

400 

769 

.458 

0 

.46123 

401 

771 

.387 

0 

.46151 

402 

773 

.315 

0 

.46171 

403 

775 

.244 

0 

.46163 

404 

777 

.172 

0 

.46167 

405 

779 

.101 

0 

.46222 

406 

781 

.029 

0 

.46290 

407 

782 

.958 

0 

.46355 

408 

784, 

.886 

0 

,46413 

409 

786, 

,815 

0, 

,46439 

410 

788, 

,743 

0, 

,46463 

411 

790. 

,672 

0, 

, 46507 

412 

792. 

.600 

0, 

.46578 

413 

794. 

.528 

0. 

46716 

414 

796. 

457 

0. 

46973 

415 

798. 

385 

0. 

47269 

416 

800. 

314 

0. 

47414 

417 

802. 

242 

0. 

47344 

418 

804. 

171 

0. 

47130 

419 

806. 

099 

0. 

46875 

420 

808. 

028 

0. 

46692 

421 
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ASTMDSP  [Fortran  version] 


4/22/88 


PAGE  1 


*****spectrum*****  modelcIO? 
thickness  is  0.0603cm 

T (1038. 7)  =  0.5033  residual  std.  dev.=  0.152E-03 

7(1107.9)  =  0.3939  residual  std.  dev.=  0.844E-04 

T(1259.8)  =  0.5110  residual  std.  dev.=  0.233E-03 

BASELINE  TRANSMITTANCE  AT  1107.9  WAVENUMBERS=  0.5057 
===>    ALPHA  (BASE)  =  0.877       ALPHA  (PEAK)  =  4.51057 

NET  ABSORBANCE  PEAK  HEIGHT  =  0.095142 

*  OXYGEN  CONTENT  IS  19.68  PPMA  * 


♦♦♦♦♦spectrum*****  modeldl2 
thickness  is  0.0588cm 

T(1040.5)  =  0.5221  residual  std.  dev.=  0.192E-03 

T(1107.8)  =  0.4015  residual  std.  dev.=  0.255E-03 

T (1259. 8)  =  0.5346  residual  std.  dev.=  0.242E-03 

BASELINE  TRANSMITTANCE  AT  1107.8  WAVENUMBERS=  0.5259 
===>    ALPHA  (BASE)  =  0.336       ALPHA  (PEAK)  =  4.33700 

NET  ABSORBANCE  PEAK  HEIGHT  =  0.102129 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦* 

*  OXYGEN  CONTENT  IS  21.99  PPMA  * 
♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦***»* 


*****spectrum*****  modeld26 
thickness  is  0.0552cm 

7(1039.3)  =  0.5151  residual  std.  dev.=  0.189E-03 

T (1107. 8)  =  0.3974  residual  std.  dev.=  0.146E-03 

T(1259.8)  =  0.5241  residual  std.  dev.=  0.241E-03 

BASELINE  TRANSMITTANCE  AT  1107.8  WAVENUMBERS=  0.5179 
==>    ALPHA  (BASE)  =  0.593       ALPHA  (PEAK)  =  4.78783 

NET  ABSORBANCE  PEAK  HEIGHT  =  0.100500 

**♦*♦♦*♦♦♦♦*♦♦*♦*********»*«*«* 

*  OXYGEN  CONTENT  IS  23.21  PPMA  * 
♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦******* 
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ASTWDSP 


4/22/88 


PAGE  2 


*****spectrum*****  modeld36 
thickness  is  0.0476cm 

T(1039.3)  =  0.5264  residual  std.  dev.=  0.172E-03 

T (1107. 9)  =  0.4355  residual  std.  dev.=  0.995E-04 

T(1259.8)  =  0.5357  residual  std.  dev.=  0.207E-03 

BASELINE  TRANSMITTANCE  AT  1107.9  WAVENUMBERS=  0.5292 
==>    ALPHA  (BASE)  =  0.303       ALPHA  (PEAK)  =  3.82890 

NET  ABSORBANCE  PEAK  HEIGHT  =  0.072947 


*  OXYGEN  CONTENT  IS  19.00  PPMA  * 


♦♦♦♦♦spectrum*****  modeld37 
thickness  is  0.0602cm 

T(1041.2)  =  0.5117  residual  std.  dev.=  0.146E-03 

T(1108.0)  =  0.4043  residual  std.  dev.=  0.168E-03 

T(1259.8)  =  0.5224  residual  std.  dev.=  0.211E-03 

BASELINE  TRANSMITTANCE  AT  1108.0  WAVENUMBERS=  0.5149 
==>    ALPHA  (BASE)  =  0.624       ALPHA  (PEAK)  =  4.12879 

NET  ABSORBANCE  PEAK  HEIGHT  =  0.091619 


*  OXYGEN  CONTENT  IS  18.87  PPMA  * 
******************************* 


*****spectrum*****  modeldSl 
thickness  is  0.0596cm 

T(1040.5)  =  0.5201  residual  std.  dev.=  0.164E-03 

T(1107.9)  =  0.4033  residual  std.  dev.=  0.102E-03 

T (1259. 8)  =  0.5301  residual  std.  dev.=  0.249E-03 

BASELINE  TRANSMITTANCE  AT  1107.9  WAVENUMBERS=  0.5232 
=>    ALPHA  (BASE)  =  0.405       ALPHA  (PEAK)  =  4.20695 

NET  ABSORBANCE  PEAK  HEIGHT  =  0.098407 

♦♦♦**♦*♦♦*♦♦♦*♦♦♦♦♦**♦♦*♦***♦♦♦ 
*  OXYGEN  CONTENT  IS  20.74  PPMA  * 
******************************* 
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ASTMDSP  [Pascal  version]  6/29/88 


PAGE  1 


sample  file  name  modeld07.dat 

thickness  (mils)  2 . 3740157480E+01 

tbase  6.0573972500E-01 

tpeak  3.9389457650E-01 

abase  8 .760629531 OE-01 

apeak  4.5103899343E+00 

anet  3 .634326981 2E+00 

uncertainty  1 .07353501 29E-02 

standard  deviation  5 . 2676600221 E+02 

ppm     (6.28)  1 .9683573442E+01 


sample  file  name  modeld12.dat 

thickness  (mils)  2.3149606299E+01 

tbase  6.2596866007E-01 

tpeak  4.0144105833E-01 

abase  3 . 3447332586E-01 

apeak  4 . 3352775874E+00 

anet  4. 000804261 5E+00 

uncertainty  1 .5548042721 E-02 

standard  deviation  5 . 2676600221 E+02 

ppm     (6.28)  2.1985050762E+01 


sample  file  name  modeld26.dat 

thickness  (mils)  2 . 1732283465E+01 

tbase  5. 1786107475E-01 

tpeak  3.9734250844E-01 

abase  5.9342275614E-01 

apeak  4.7850185869E+00 

anet  4. 1915958308E+00 

uncertainty  1 . 2867481 063E-02 

standard  deviation  5 . 2676600221 E+02 

ppm     (6.28)  2.3183221817E+01 
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